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We report on a study of a finite system of classical confined 
particles in two-dimensions in the presence of a uniform mag- 
netic field and interacting via a two-body repulsive potential. 
We develop a simple analytical method of analysis to obtain 
ground state energies and configurations. We prove analyti- 
cally the minimum energy configurations are independent of 
the nature of two body interaction and the magnetic field. In 
particular we prove that the first transition from a single shell 
occurs when the number of particles changes from five to six. 
These results are exact. 

PACS numbers: 02.60. Cb, 71.10.-w 



There has been considerable progress in recent years 
in the study of electrons in quasi-two-dimensional sys- 
tems both experimentally and theoretically There 
are several examples of these systems but the most im- 
portant one from our point of view is that of electrons 
in a quantum dot. There have been several studies on 
the ordering and transitions of charged particles in two- 
dimensions ||-^. In a recent Monte Carlo study, Be- 
danov and Peeters ||](see also Bolten and Rossler 
have analysed the classical ground state of a system of 
confined charged particles interacting through Coulomb 
interaction. By minimising the classical energy, they ob- 
tain numerically the shell structure in a cluster of N- 
particles. They have systematically listed the shell struc- 
ture in a "Mendeleev" table for < 52, and for a few 
large clusters. Similar results are available also for loga- 
rithmic two-body interaction 

In this paper we consider such a cluster in which the re- 
pulsive two-body potential is either a power-law( or loga- 
rithmic). We devise a simple analytical method to obtain 
the classical ground state energy in two steps. First we 
minimise the energy for a fixed total angular momentum 
(which is conserved), J and then minimise this energy 
with respect to J. This has two advantages. It reflects 
the quantum degeneracy of the lowest Landau level for 
electrons in a uniform magnetic field in the absence of 
any interaction, at the classical level itself and secondly 
it allows one to do the second step minimisation over the 
quantised values of the angular momentum. Here how- 
ever we restrict ourselves to classical analysis only and 
derive some exact results analytically. 

In particular we show that: a) the configurations min- 
imising the energy are independent of the parameters of 
the Hamiltonian (eg magnetic field) so long as the re- 
pulsive two body potential falls off as a power-law (or 



vary logarithmically) with the relative distance. Only 
the overall length scale is sensitive to these details. 2) 
Two special configurations in which all the N particles 
are on a circle (referred to as 0) ^-^d the one in which N 
- 1 are on a circle with one at the center (referred to as 
Q ) are always (local) minimum energy configurations. 
We give exact analytical expressions for the correspond- 
ing minimum energy for all N . The Q ^i^s lower energy 
for iV < 5 while has lower energy for iV > 6. This 
geometric transition is the first one to occur and is inde- 
pendent of the precise form of the repulsive interaction. 
3) While it is known numerically that for iV > 9 and 
for Coulomb potential |^,^ the minimum energy config- 
urations exhibit approximate multi shell structure, the 
special configurations provide an upper bound on the 
minimum energy for a whole class of interactions that 
we consider here. 

The classical system we are interested consists of N 
particles confined in an oscillator potential in a uniform 
magnetic field and interacting via a two body interaction. 
The Hamiltonian of such a system of particles given by. 
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where and pi denote the position and momentum vec- 
tors of the i-th particle. The vector potential, for a uni- 
form magnetic field is given by. 
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where ojl is the Larmor frequency and rij = ri — fj . The 
power v (positive) is kept arbitrary. In what follows we 
also comment on case when the two-body interaction is 
of the form —(3 J^i j(-ii) ^og{r^j/p^) which is repulsive for 
rfj < p^. This may be closer to the real situation as 
also the i/ = 1/2 (Coulomb) in the case of electrons in a 
quantum dot [Q. Recently, the v = 1 case has also been 
analysed in detail The Hamiltonian can be written 
in terms of dimensionless units by introducing a length 
scale I = (hj (muj) which is the basic oscillator length. 
All distances are measured in terms of this basic length 
unit. Note that the h is introduced only as a convenience 
so that the energy is measured in units of hto and does 
not have any other significance as in the quantum case. 
The momenta are measured in units oiTi/l. 

The new Hamiltonian in these scaled units, but keeping 
the same notation, may be written as 
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where ji — r,; x p,;, a = 77 and g — ■^(O^'^- Unless 
otherwise mentioned the summations run from 1 to 
hereafter. While the original coupling constant (3 was 
dimensional the new coupling constant g is dimension- 
less. Hereafter we assume all the energies are measured 
in units of fiui and do not write the units explicitly. 

For the first step of the minimisation we introduce the 
function, 



F = H + \{J -Y^j,), 
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where ji are the single particle angular momenta and A 
is the Lagrange multiplier which enforces the constraint 
J = Y^iJi- Setting 5F — 0, gives the necessary equations 
to determine the equilibrium configuration in the phase 
space, 
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These are the basic set of equations. Any solution to this 
set of equations describes an equilibrium configuration 
but not necessarily the one with the minimum energy 
for the given J. The case of logarithmic interaction is 
obtained by simply setting the power 1^ = and by set- 
ting the prefactor to 4g instead of Agv in eq.(|7|). With 
this proviso, unless otherwise mentioned, all subsequent 
equations also reproduce the log case. 

First we present a qualitative but a general analysis of 
these basic set of equations. To make the analysis simple, 
we introduce an auxiliary variable, 
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The total angular momentum may now be written in 
terms of this auxiliary variable as 



J = 2_^fi X Pi = -{a + X)<j), 
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where we have made use of eqs.(|[||). It is convenient 
to express r,; = Rsi with R being a common scale fac- 
tor which may be taken to be the radius of the farthest 
particle, say the A^*'* one. Therefore 
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since sf^ = 1. Using eq.(0) and eliminating A dependence 
using eq.(H) we have, 
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Taking scalar product with Si and dividing both sides by 
si (7^ 0), we get. 
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Note that the LHS is independent of the particle index 
i. Thus we have iV — 1 independent of equations of the 
type 
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Further by taking the cross product with Si and dividing 

by Si, we get. 
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which provides a further set of N conditions on the inter- 
nal coordinates Si . Notice that these conditions are man- 
ifestly scale invariant. Together eqs.( p^ , p^ provide the 
2N — 1 necessary equations for determining the s^ and 
the angles Qi. Notice that these determining equations 
are completely independent of a, J and g. We have there- 
fore the result that (si, S2, sat-i, 0i, ^2, ^Af) are in- 
dependent of the magnetic field (a), the total angular 
momentum J and the interaction strength g. These pa- 
rameters, however, determine the overall scale R through 
eq.(|l^). In fact since s^ = 1, the corresponding equa- 
tion may be taken to be the determining equation for R 
in terms of the parameters of the Hamiltonian, 
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Note that A{i', N) introduced above is a function of 1/ 
and N only. Thus we have a very general result that the 
geometry or the shell structure of the equilibrium con- 
figuration is independent of the parameters of the Hamil- 
tonian which only restrict the overall size of the system. 
The shell structure, however, depends on the nature of 
repulsive interaction through the parameter i/ but not its 
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strength. The above analysis is vahd even if any one of 
the Si = 0, that is one particle being at the origin of the 
coordinate system in which case we have two equations 
less ( not more than one particle can be at the origin). 

The energy of the equilibrium configuration can be eas- 
ily computed by noting that the auxiliary variable (f) de- 
fined in eq. (^) is related to the two body potential energy 
by, 
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where the RHS is proportional to the potential energy 
due to interaction. Since the RHS and are positive 
definite we have the condition (1 + > J^. Inter- 

estingly taking further gradients of the above equation 
again yield conditions on </) which violate the above in- 
equality. Thus the minima described by the equations 
®il 0) must be isolated. We have for the energy at the 
extrema, 
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where 4> = B?'4> and </> is independent of J. This then is 
the energy of the equilibrium configuration in a given J 
sector. For the logarithmic case, the first line of eq.(p^) 
has the second term corresponding to the logarithmic 
interaction and the second line is not valid. 

In order to find the global minimum we now minimise 
the energy with respect to J and set dE/dJ = 0, that is, 
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Differentiating the R w.r.t. J in eq.(|T^) yields dR/dJ 
which when substituted in the above equation gives the 
J value and the corresponding global minimum energy, 



absence of the magnetic field as it should be. The ex- 
pressions given above, though, are valid independent of 
the geometry of the clusters and are exact ( for approxi- 
mate solutions see eqs.(8,9) in ref. [Q for the special case 
of Coulomb interaction) . 

The geometry of the clusters or shells are dependent on 
(j) and A(v, N) which are as yet unspecified. In general 
the equations for the equilibrium configurations admit 
many solutions (which are isolated as remarked earlier) 
for a given N and i^. In particular there are two special 
configurations which are always solutions viz (i) all the 
A'' particles are on a circle , Q ^^^d (ii) iV — 1 particles 
are on the circle with one particle at the center, Q . For 
these two cases only the over all scale factor R is to be 
determined. The angles, 0ij/2, are simply multiples of 
tt/N and tt/{N — 1) respectively. These however need 
not be minimum energy configurations for a given N and 
i/. In fact it has been numerically proved that for < 5 
the circle configuration is indeed the minimum energy 
configuration where as for 6 < iV < 8 it is the circle-dot 
which is the minimum energy configuration in the case 
of Coulomb interaction ( i/ = 1/2 ). Multiple shells start 
forming for N > 9. In what follows we prove analytically 
that the first transition which occurs for N from 5 to 6 
is independent of v. The case of circle and circle-dot is 
particularly simple since there is only one scale involved. 
That is all s-f — l,i — 2, ...,N and sf — I for the circle 
and sf — for the circle-dot. 

For the circle case, we have, 
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and therefore the energy is given by, 
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In the case of circle-dot, we have. 
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where 
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and A{v^N) is defined before in eq.(|l5|). An important 
point to note here is that the minimum energy E is in- 
dependent of the magnetic field and its dependence on 
g is explicit. The dependence on N and v is however 
involved. The angular momentum J at minimum of en- 
ergy depends on the magnetic field and is zero in the 
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since there are now — 1 particles on the circle and 
therefore the energy is given by. 
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The extra 1 on the RHS is due to the contribution of the 
particle at the center. To ascertain which of these two 



3 



configurations ^^'^ O has lower energy it is sufficient 
to look at the ratio, 



where, 



N 



N -1' 



X 



N 



22^4 



(27) 



N-1 



1 



^1 sin2''(f ) 



(28) 



Note that in general the ratio / depends only on N and 
V. Obviously the circle is a lower energy configuration iff 
/ < 1. We claim that, for all u > 



f{u,N) < 1 for < 5; 
f {v, N) > 1 for N > 6. 



(29) 



Further the function /(i/, N) crosses unity exactly once 
for N between 5 and 6 and nowhere else. This result 
can be easily seen for = 1 since in this case A^'' = 
{N^ Therefore, 
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which reproduces the claims made above, for ly = 1. The 
general proof of these claims for all v is some what in- 
volved and will be published elsewhere. We sketch the 
arguments for the case v «\^ where one may use the 
expansion a"^ ~ 1 + ;/log(a). Using this Aa? may be writ- 
ten as, 



out that the difference is precisely given by g[iN- There- 
fore the first geometric transition also occurs for the log 
case exactly as in the power-law case. 

To summarize, we have proved in general that the or- 
ganisation of many body clusters in two dimensions into 
shells is a robust phenomenon independent of the nature 
of the repulsive two-body interaction and also indepen- 
dent of the Hamiltonian parameters but dependent only 
on the number of particles in the cluster. In particu- 
lar we have analytically proved that the first geometric 
transition for the ground state from circle to circle-dot 
configuration occurs after iV = 5. The robustness of this 
transition seems to emerge purely from the number the- 
oretic properties of the ratios of the energies (difference 
in the log case) in these two configurations. It is an open 
question if this is due to some hidden symmetry proper- 
ties. We have also done numerical simulations for larger 
N and for various values of v. We find that the shell 
structure found by Bedanov and Peeters |^ in the case 
of Coulomb interaction is valid for all v in general except 
that for larger v the shell description is valid only ap- 
proximately. The equilibrium configurations correspond 
to isolated minima. This fact should be useful in calculat- 
ing quantum fluctuations about the minima. The details 
of these investigations will be published elsewhere. 
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where 
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where we have used the identity O^i^ sin(^) — 
Substituting for Xn in /, we have 
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N{N - 1)' 
where 

fiN - [N{N - 3) log(iV) - (iV - 1){N - 2) log(iV - 1)]. 

(34) 

is independent of i/. It is now easy to see that /iat is neg- 
ative for N < 5 and positive otherwise. Hence the claim. 
In the case of logarithmic interaction, the transition can 
be seen more easily by taking the difference of energies 
for circle and circle-dot since this difference is indepen- 
dent of the arbitrary scale p in the interaction. It turns 
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